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A ataody  atat*  nodal  (or  (luortnr  aton  rrconbinatlon  for  tan^raturaa 
laar  than  900  K on  a natal  catalytic  turlaca  uaa  davtlopad  baoa4  on  tha 
•Ideal  hatarogaoaoua  raaatton  nacHanlan.  TKa  nodal  la  tou|»rlaa4  of  thraa 
dlatlnct  procvfiaaa:  (1)  tha  bomijt\ft  yrocaa*  of  gaa-phata  fluorlna  atoaw 

to  tha  natal  aurfaca,  (7)  tha  collision  and  auhscauont  raconhtnatfon 
yrocaaa  of  gaa-f»hasa  aton?'  and  aur I acr-hondad  flucn  lna  atona  and  (})  tha 
rata  yrucwsa  at  which  the-  bonding  end  rrfoPHn»», i<?^  taVa  fon''lda"-- 
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ant*  CkAMiriCAt 


lac  the  Kurfac*  coverage.  The  ar^el  waa  eeegar»t!l  with  aclwal  fluorine 
reetwliinatien  gate  taken  for  « iiickel  tuhe  In  which  a gee  nlRturn  of 
atonic  fluorine,  nolrrular  fluorine,  and  argon  flovtHif.  The  reeulta  of 
thla  ceaiporlnoa  vallt^ateg  the  npgel  and  eatended  ihn  vartficacloa  of  the 
RI4m1  MchanlM  for  heteregewona  racodhttiattMi, 

Clef  Mtat tana  watng  thn  wall  reeoalitiMU ion  icmNl  were  Wkle  Inr  a 
t«N»-dli«Miatomal  all*  nogale  of  the  t^gd  i^ieed  la  flnorlae  ahcnlcal  laaera, 
natag  a t«o<gliHM»atonal  "'hMiirall|r  reacting  leadwar  houndary  layer  cea- 
gater  coda,  iht  reewita  t t.heae  cregutatlotia  ahowtrd  that  the  wall  recei 
hlaatle*  auhatantlally  afi acted  the  tungeratare  grefllce  ee  well  e«  welt 
heat  iranefar  la  cbealcal  laaet  iK>taiea<.  The  aMuatt  nf  atiailc  fluorine 
degilaled  through  wall  racoedilnatlon  wa^  ihown  to  ha  on  the  order  nf  a 
third  of  the  laltlcl  flnerina  concwntratlcm  lor  the  naaala  caaf Igaratlen 
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Tht  orl|[in  of  this  diosortotion  dotot  back  to  a trip 
t wmdt  to  tka  Ifoapaat  taberatory  ia  March  IfTS.  At  that 
tlaa  tha  problta  of  flooriao  rocoabination  toosad  to  m to 
ho  am  oalroMOly  praaaiiig  oo'^»  and  ooa  ohich  1 thought 
uould  cortaialy  bo  aolvod  by  tha  tioa  I «rai  raady  to  bagla 
•y  rotaarch,  I rotumod  to  AFtt  aod  bogaa  to  proparo  oy* 
•otf  for  oubarklag  on  a parauotric  if^y  of  boat  transfer 
of  foots  to  fluorlno  rocaibinatlon  on  tha  walls  of  tha 

••all  nottlts  in  a fluorlno  choolcal  laser  usini  soiaoono 
olao's  rocoibination  nodal  (which  t assunod  would  bo  in 
aalatonco  when  1 naodad  it). 

Vhlla  t was  back  at  AFIT,  an  lbtara»ting  thing  took 
ptoca  In  fluorLst  ehtnical  lastr  circles.  It  was  tha  di«* 
cotrary  of  tho  infonous  *^t  Chain*'  fluorina  laser  reaction. 
In  effect  tho  Hot  Chain  was  adwartisod  as  btlng  sbl#  to 
produce  excited  HF  and  OF  using  a reaction  which  required 
only  snail  anounts  o^  atonic  fluorine  and  large  anounts  of 
noltcular  fluorine.  The  effect  of  this  "discovery"  was  to 
put  an  and  to  the  urgen*:y  of  developing  a fluorine  wall 
reconblnat ion  nodel.  Thus,  by  the  tine  I arrived  at  the 
Neepons  Laboratory  to  bagln  aiy  paranetric  study,  nothing 
had  been  done  on  developing  a reconblnat  Ion  nodel.  Tht 
■odet  was  still  considered  an  Interesting  problen  bu*  no 
longer  an  urgent  one.  So,  I began  the  work  of  developing 
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tiM  rtcpiAinatlon  aodal  p?«8aat«4  htrt.  In  tka  aaiM  tlaa 
fraaa  s<mm  of  «arly  toats  uilng  tha  Hot  Chain  warn 
%•!•§  taHtortakan  with  nagatlTo  raaulta  bnt  high  kapaa. 

Akaat  ala  aaatka  ago  paapla  bagan  raalita  that  tba 
Nat  Chain  waa  not  all  it  ka4  baaa  aatcwad  H.  ba«  S«44aal]rt 
atll  racoabl nation  araa  again  an  inpoTtaot  and  platting 
pfoblan  and  t bad  paapla  aibing  for  ay  nodal ^ In  tkit 
ragardt  I found  tba  work  naat  ancitingi.  aapacially  tinea, 
bacauta  of  tka  tint  Involvad  in  parfaming  tba  raaaarck, 
naat  diaaartation  topica  aia  not  of  tuck  lanadiata 
iaponanct.  So.  at  I bagin  tbanking  tkoaa  indioiduala 
ubo  aidad  no  in  nv  work,  t would  firat  lika  to  tkank  who* 
avar  it  waa  who  firat  auggaatad  that  tka  antwar  to  fluorina 
wall  raconb) nation  waa  tka  Hot  Chain. 

In  all  aarieuantaa,  t would  lika  to  tkank  all  tkota 
at  tka  Vaapona  Laboratory  %rko  wara  of  auck  kalp  to  na,  in 
particular,  LtCol  Jay  koland,  LtCol  David  Ericaon,  Ma) 

Dick  Harris , a apodal  thanks  to  LtCol  Dick  kobards  for 
his  patiant  «ar,  and,  cartainly  not  loast,  Mrs  Eatkar 
Castor  for  typing  this  diasortatlc  i.  1 an  also  indabtod 
to  savaral  nonbors  of  tha  APIT  faculty,  in  particular,  ny 
advisor  Dr  Janos  Hitchcock  for  Ms  halpful  suggaationi  and 
connonts  and  Dr  Andraw  Shiaa  for  hit  confidanca  in  na. 
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Abitract 

A «t*ady  stat*  for  floorino  ttn*  rocooblnation 

for  ttaptraturos  lt«s  than  109  C on  a natal  catalytic  aur- 
faca  vat  davalopad  batad  on  tha  Rldaal  htiarogtnaout  raac* 
lion  nachanita,  Tht  nodal  It  coaprltad  of  thraa  dlatinct 
procattat:  (1)  tha  bonding  procatt  of  gat*pha«a  fluorlna 

atont  to  tha  natal  turfaca.  (2)  tha  collitlon  anj  <.*’bta* 
^uant  raconblnat ion  proce^t  of  gat-phata  atom  end  turfaca- 
bondad  fluorine  atom  and  (5)  tha  rata  procatt  at  vhlch 
tha  bonding  and  raconbination  taka  place  contlderlng  the 
turface  coverage.  The  nodal  vav  conpared  with  actual 
fluorine  raconbination  data  taken  for  a nickel  tube  in 
vhlch  a gat  nixture  of  atonic  fluorine,  nolecular  fluo- 
rlna, and  argon  flovad.  Tha  ratults  of  this  conparlton 
validated  the  nodal  and  extended  th*  verification  of  the 
Ridaal  nechanisn  for  beterogeneous  raconbination. 

Important  characteristics  of  fluorine  vail  raconbina- 
tion were  denonstrated  through  use  of  the  nodel.  The 
raconbination  coefficient,  defined  as  the  ratio  of  the 
raconbination  resulting  vail  collisions  to  the  total  vail 
collisions,  vas  found  to  be  a strong  function  of  tenpera- 
ture.  The  recombination  coefficient  dropped  rapidly  from 
a naxinum  of  approxinaiely  10  at  400  K to  approximately 
% X 10*  at  700  I due  to  thermal  desorption.  The 

alv 


rtco«bin«tLon  co«ffici«n»  wat  also  found  to  It  stron|iy 
deptndtnt  on  tht  diluent  of  the  gas  nixture  for  lower  tea* 
paraturts.  An  argon  diluent  was  found  to  greatly  Halt 
tint  recoaibinat Lon  coefficient  at  teaperatures  below  4C0  K. 
Neliua  is  predicted  to  have  auch  less  of  a llaltiag 
influence  when  used  as  a diluent. 

Coaputatlons  using  the  wall  recoabinatlon  aodtl  were 
•ade  for  a two*diaensional  slit  nottlt  of  the  type  used  In 
fluorine  rhenlcal  lasers,  using  a two*dlaenslonal  chen;* 
tally  reacting  laninar  boundary  layer  computer  code.  The 
results  of  these  coaputations  showed  that  the  wall  recoa* 
bination  substantially  effected  the  teapersture  profiles 
as  well  as  wall  heat  transfer  In  chenical  laser  nodes, 
The  aaount  of  atonic  fluorine  depleted  through  wall 
recoabinatlon  was  shown  to  be  on  the  order  of  a third  of 
the  initial  fluorine  concentration  for  tht  noztle  con* 
figuration  run. 


A »«ODEL  FOR  FLUOR?NE  i-ECONBINATION 
AT  A METAL  SURFACE 

t . Introduction 

incItTound 

For  «o»«  r«nr»  nou  it  hn»  b«<>n  thut  cheairal 

Insert  hairn  a nuch  Muhar  perfontance  potential  than  gat 
dfnanlc  laser*,  (IDL,  in  teras  of  voluaetric  efficiency, 
cHtnical  efficiency,  and  gain  (Ref  1).  Recause  of  the 
siupler  technology  involved,  GDI’s  were  the  first  lasers 
to  he  esplored  at  the  development  level  for  the  production 
of  a high  power  laser.  Research,  however,  has  continued 
on  cheaical  lasers  tc  lay  the  grounti  work  for  the  second 
generation  of  high  power  lasers.  Since  19?1  considerable 
effort  has  been  expended  in  the  production  of  comprehensive 
gesdynaisic  computer  models  for  the  study  and  design  of 
chemical  lasers  (R«*f  ?)• 

Computer  modeling  of  thr  norzle  banks,  lasing  ravi- 
tits  and  diffusor  sections  is  Just  beginning  to  gain 
acceptance.  Usable  results  are  beginning  to  be  produced 
toward  the  final  goal  of  merging  the  separate  Models  of 
each  of  the  compone.nts  Into  an  overall  design  code.  Of 
major  importance  to  these  codes  are  the  phystuchemical 
inputs.  These  inputs  are  presently  being  looked  at  in 
great  detail  in  an  attrnpt  to  further  refine  the  semi- 
tmpirical  models, 
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One  area  which  has  been  in  question  for  a nuaiber  of 
fears  is  that  of  the  fluurir.e  wall  recombination  reaction 


in  hydrogen  fluorine  and  deuterium- fluorine  lasers.  To 
•id  to  understanding  the  effect  of  the  fluorine  wall  recon 
blnailon  on  the  ^p«ration  of  a fluorine  chemical  laser,  r'' 
Mke  up  and  operstion  of  such  a system  is  «»xanlned. 

Figure  I is  a Khematlc  representation  oi  the  general 
class  of  fluorine  chemical  laser  systems  of  interest.  The 
primary  flow  contains  the  fluorine  atom  source  to  be  used 
in  the  chemical  reaction  which  will  produce  the  excited 
species  for  lasing.  This  primary  flow  is  usually  the 
result  of  some  earlier  chemical  combustion  reaction  such 
at 


bHj,  ♦ bFj  • «►  di  • eMf  * ♦ cHe  (1 

where  the  Hj  is  pirseni  in  ?-.ifficient  quantity  to  cause  a 
large  enough  reaction  with  the  fluorine  to  raise  the  ten' 
perature  of  the  gas  to  a temperature  level  compatible  with 
the  Ptslred  no::ir  input  quantities  of  atomic  fluorine, 
which  is  for  the  most  part  present  from  thermal  dissorit* 
tion  of  the  molerular  fluorine.  The  helium  is  representa- 
tlve  of  a diluent  added  to  bring  the  properties  of  the  gas 
to  those  compatible  with  the  total  laser  flow  system, 
nottle  bank,  lasjrig  cavity  and  diffuser  section. 


Ch«fflical  lastr  Cp«r«tlo» 


Rtferring  to  Figure  2,  the  priaer^r  flow  iii  then 
•ccelereted  through  the  priaary  noisle?  of  the  notzle  bank 
reaching  velocities  ranging  fro*  Math  2 to  Mach  S depending 
on  the  particular  laser  configuration.  It  is  on  the  walls 
of  these  not  ties  that  the  fluorine  re«.o*binat  ion  tan  take 
place.  Upon  exit  fro*  the  nottles  the  prlaary  notzle  flow 
is  joined  in  the  nixing  and  lasing  cavity  by  the  adjacent 
Secondary  noitle  flow^  vonsisting  oi  either  or  Dj.  The 
jets  nix  In  the  nixing  zones  indicated  in  Figure  ! where 
the  excited  HF(v)  lasing  states  are  forned  in  reactions  of 
the  type 


• F -MF(v)  • H f2) 

The  nirrors  for*  the  optical  isvitv  and  power  is  extracted 
fro*  the  '•avSty  through  ti»e  wirror  havjnj.  psriial  trans- 
nittance. 

There  are  nany  factors  hearing  on  the  natching  of 
theoretiial  gam  data  vtth  oh^erve.!  r.speriisefit«l  gain  data 
(Ref  J).  To  nention  m few  areas  where  inadequate  woJrling 
nay  be  resportsiHe  for  the  niswatrhmg  the  following  are 
given 

a)  P'jwping  rr  u t mn  rates  for  reactions  of  the  type 

F ♦ m (v)  ♦ H (.1) 
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b)  Stiaulated  enisston  relMXAtxon  rat««  for 
reactions  o!  the  tyye 


Hf{t)  • hv  -HF(v-l)  • 2hv 


c)  Vibrat tonal ■ rotat ional  and  vibrational' 

translational  deactivation  rates  for  teactlona 
of  the  type 


HF(v)  » M •HF(v-6v)  • M 


d)  Sinfle  and  nul  t i *^uanluM>  vihr a t iona  1 > v ibra  t iona  I 
transfer  rates  for  reactions  of  the  type 


HF(v)  * HFtv')  ••Hr(v*ftv) 

• Hf(v‘  Av) 


e)  Turbulent  nising  and  diffusion  proce<se'‘ 

f)  Free  itreae  floor  ine  rcf  onbinat  inn  rate'  <‘or 
reactions  of  the  type 


2F  * M *-F^  * M 


|)  Kail  fluorine  reconbinat  io.'  rates 
It  is  unclear  at  this  tine  exactly  how  each  of  the  factors 
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directly  influence  laser  gain,  as  each  plays  an  interrelat- 
ed role  with  the  others.  Of  particular  interest  in  this 
itudy,  however,  is  item  g above.  It  is  clear  froa  Ec|ua*  iv.i 
C2)  that  the  vail  t evoab inat ion  reaction  is  a detriaentwl 
one.  The  less  r available  to  the  excitation  reaction,  the 
lest  excited  species,  and  thus,  the  less  output  j over 
possible  fron  the  laser.  Further  understanding  of  the 
posslMe  effect  of  this  reaction  has  been  pointed  out  by 
Ferrell,  Kendall,  and  Tong  (Ref  4)  and  Miilatanan,  Kurtius, 
Me  Dana  1 , and  Thorenes  (Ref  S),  Ferre  1 ei.al.  showed  that 
one  night  expect  a large  reduction  in  atonic  fljorine  and 
a pronoun  ed  gradient  in  the  atonic  and  molecular  c^vneen 
tration  profiles  at  the  e.sit  plane  of  the  primary  norrles 
in  a study  in  vhich  every  fl  torine  a too  colliding  with  the 
wall  was  assuned  to  result  in  a recombination,  Mikaiarian 
et.al.  showed  the  effect  one  niaht  expect  froiD  such  a con 
centraiioh  profile  on  the  laser  output.  In  an  attempt  to 
explain  •nonalit*''  between  calculated  and  obsersed  perform- 
ance curve'',  thev  input  a crudely  e*;tiRated  v oni  ent  rat  ion 
profile  into  a cavttv  and  la^jnn  code.  The  result 

was  that  a notMini  fert»  initial  fluorinr  i onv  ent  rat  inn  pro 
file  in  the  laser  cavity  entrance  caused  twc>  important 
effBCts.  reduction  in  the  gain  due  to  a decrea'.e  of 
reactants  and  axial  shifting  of  the  excitation  reactions 
down  stream  fron  the  exit  plane  of  the  norrlrs  (ree  Figure 
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5).  Thus  A - iind  understaml Srvg  if  <his  react  iojv  i of 
inportence  in  any  attenpt  tcf  •aoi.lel  fluorine  i heaical 

lAfem . 

H*t«rOi|!eneQu»  Re yoiwhinat  ton . In  general  the  litera 
turtr  containit  abundance  nf  vork  dealing  with  the  prob 
It«  of  letonHnaiion  at  a catalvtic  mi?  face.  For  the  «o»t 
part  the  work  is  exper tv^ental . Varicus  «echani**s  for  the 
Mall  ref owbinat ion  process  have  been  proposed  in  an  attenpt 
to  correlate  the  data.  ■'U  these  aechanisns,  three  appear 
to  he  possible,  tvo  involving  onlv  surface  adhered  atows 
proposed  bs'  HinsheNood  and  one  involving  a reaction 
between  A gas  phase  and  a s?irfare  adhered  ato«»  bv  Rideal 
(Kef  b,*)  These  ne<han(sms  are  brien>  described  below 

(1)  In  the  Minshe  laood  rve » *:hbo)  J nj;  adhered  ator  aoJel 
it  is  proposed  that  the  n«*cbanis»  for  r et  owl  i nut  ion  i<  the 
adsorption  of  the  reactants  at  neighboring  surface  sites, 
Kecaiise  of  the  nevessitv  of  neiph?*oring  collisions,  this 
Hodel  predSi  ts  .?  set  ond  orvier  re;uTion  rate  vifh  lesprct  to 
pressure  at  all  te^'perat  «u  os  . 

(*’)  In  the  Mifshelvv'oi)  'rigrati  -n  of  adhered  ator 
Node)  it  is  propo'se.i  that  the  mechanisn  foi  re;  onV  |nat  ion 
is  the  slov.  adsorption  of  :e.ait.irt  sjM;ir  amwhere  on  the 
surface  am!  rapid  surface  ^igtaliorn  of  the  adsorbed  reac 
tants  folloved  bv  recoi^binat  ion . Since  ever;'  collision 
which  prodmes  an  adhered  ator*  \ .;»n  viiuse  » reaction,  this 
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C.  CalcuUt  ionx  Maue  with  Species  Profiles 


Shown  Above 


nod#!  predicts  a first  order  reaction  rate  with  respect 
to  pressure  at  aHI  teaperatures . 

())  In  the  Rideal  aodel  It  Is  proposed  that  the 
•echanlsM  for  recoabination  Is  a reaction  by  collision  uf 
e gas  phase  atoa  with  an  adhered  surface  atoa.  This 
aechanisa,  unlike  the  two  previous  ones,  predicts  both  a 
first  and  second  order  reaction  with  respect  to  pressure 
over  a range  of  teaperature. 

The  currently  favored  aechanlsa  is  the  Rideal  aodel 
since  the  expertaental ly  observed  pressure  dependence  is 
first  order  at  low  teaperatures  changing  to  second  order 
at  higher  teaperstures  for  slow  wall  recoabinat ion  reac- 
tions. This  is  the  saae  teaperature  dependence  predicted 
by  the  Rideal  aodel  (Kef  6,R). 

With  the  use  of  the  Rideal  aodel  and  the  assuaption 
of  an  Airhcnius  rate  process,  the  various  catalytic  wall 
reconbinat ion  reactions  have  been  quantified  using  the 
eaperinental  d«ta  for  that  reaction.  Notably,  howerer,  a 
theoretical  ba«iis  for  explaining  soar  experimental 
anomalies  is  very  veak  (Ref  9),  Because  of  the  lack  of 
exper  iirent  al  data,  catalytic  wall  reconbinat  ion  models  for 
fluorine,  and  in  fact  .ill  the  Halogens,  do  not  exist. 

Review  Available  Data  for  Halogen  Ralj^  Recewbinat  ton. 
Very  little  data  on  Halogen  surface  recoahinat ion  were 
found  in  tbe  open  literature.  Work  done  by  K.  A.  Ogryilo 
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(R«f  10,  11)  indicated  that  approxinately  one  in  ten  wall 
colH«ionik  resulted  in  a recombination  for  chlorine  and 
bromine  atoms  with  a metal  surface  at  room  temperature. 
Ofryilo  reported  that  fluorine  was  too  difficult  to 
handle  and  thus  no  fluorine  rates  were  measured. 

V.  Valance,  B.  Birang,  and  D.  Maclean  (Ref  12) 
reported  that  approximately  one  in  10^  wall  collisions 
resulted  in  a recombination  for  fluorine  on  quartz  at 
room  temperature.  No  work  was  completed  by  this  team  on 
interactions  of  fluorine  with  metal  surfaces. 

A carefully  prepared  set  of  unpublished  data,  however, 
was  found  to  exi^t  for  fluorine  reconh inat ion  on  a nickel 
surface  and  was  nadr  available  for  this  studv  by  Dr. 

Casper  J.  Ultee  of  United  Aircraft  Research  Laboratories 
(Ref.  15). 

Object  t vc«-  and  of  the  rresent  ‘ udy 

This  study  had  three  objectives;  (1)  Develop  for  the 
first  time  a qu.snt  ii  tat  i vc  model  for  fluorine  recenhinat  ion 
at  three  rrt*l  surfaces  vooronly  used  in  laser  nortle  con 
struction  • copper,  Rold,  and  nickel.  (2)  ("ompare  the 
model  with  existinR  fluorine  wall  reconbinat  ion  data. 

(5)  Use  the  model  to  e.vas»ii)e  the  effect  of  "real*'  fluoiine 
wall  reconbinat ion  (nudel  verified  by  experimental  data)  on 
the  flow  in  a typical  chcmisai  laser  no;  :1c. 
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Xhm  d«v«lopa«nt  of  th«  fluorine  well  reciHiblnation 
•odel  Is  presented  In  Chapter  II.  The  approach  was  to 
adapt  the  Rldeal  aechanita  as  the  doalnant  reaction  •ech- 
aalsM  (consistent  with  known  recowblnat ion  principles), 
dacoapose  this  aechanisa  into  its  coaponent  physical 
processes,  and  quantify  each  of  these  processes  froa  basic 
principles.  This  differs  froa  the  conventional  approach  of 
proposing  a aechanist  and  adjusting  its  bulk  paraacters  by 
choosing  the  best  fit  to  experiaental  data,  by  taking  the 
basic  principle  approach,  the  aodel  is  able  to  take  into 
account  the  detiil«d  environaent  in  which  the  reaction 
takes  place.  This  leads  to  the  inclusion  of  diluent 
effects,  which  have  not  o*en  carafully  treated  n wall 
reconbinat ion  theories  for  any  gas. 

Chapters  III  and  IV  deal  with  objective  two,  the  coa> 
parison  of  the  nodel  with  the  Ultee  reconbinat ion  data. 
Chapter  111  develops  the  analytic  tools  necessary  for  this 
conpariscn  and  the  results  are  presented  in  Chapter  IV. 

Chapter  V deals  with  objective  three,  detereining  the 
effect  of  "real’'  fluorine  wall  reconbin.it  ion  o.;  the  flow  in 
a particular  chenical  laser  norrle.  A description  of  a 
steady  state,  two-dinen$ional  chenirally  reacting  laninar 
boundary  layer  conputer  code  developed  for  this  study  is 
included.  Finally,  the  results  of  the  nozzle  analysis  are 
presented. 
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II.  Fluorine  >fall  Recotablnat  ion  Model 


Til*  lUdeal  Mchanian  waa  adapted  aa  the  doninant  raac 
tiofi  Mchanlaa  becauat  of  pravioua  verification  for  slov 
reaction  reilaea.  The  tem  '‘aodier*  la  taken  here  to  naan 
the  quantification  of  the  terns  in  the  rate  equations 
daacribini  *hc  proceasei  of  the  wall  reconbination  nacha* 
nlan.  The  nodel  seeks  to  adequately  describe  a reactive 
interaction  between  a g»%  phase  fluorine  aton  and  a surface 
adhered  fluorine  aton  to  produce  a fluorine  gas  molecule 
in  un  attempt  to  extend  verification  cf  the  Rideal  necha- 
nisn  to  reaction  reiimes  which  differ  from  those  pre- 
viously verified.  Gas  kinetics  were  examined  to  describe 
the  gas  phase  aton  behavior  at  the  surface.  The  aton- 
wall  intera':tion  phenomenon  was  examined  to  establish  the 
surface  aton  concentration.  And  finally,  the  rate  equa- 
tions were  written  and  solved  for  the  steady  fiate  case  to 
fully  establish  the  reaction  model.  A summary  of  each  of 
these  effort«i  is  given  in  the  fol loving. 

Gas  Phase 

Kinetic  theory  was  used  to  derive  the  rate  of 
inpingenent  of  atoms  on  a surface  and  probability  distri 
billion  for  the  normal  component  of  the  velocity  of  the 
impinging  atoms.  Th*  differential  flux  v>f  atoms  on  a 
surface  for  a bolttmann  distributed  gas  is  given  by 
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(Kef  14) 


dA 


dVjdV^dVj 


(•) 


where  n Is  the  nu«b«r  density  of  etons,  Vp  and  are 
the  cartesian  velocity  con|»onents  of  the  atoes.  where 
if  the  nornal  coMponent , ■ is  the  ness  of  the  aton,  k Is 
the  Boltznann  constant,  v>nd  T Is  the  tenper 4ture ..  The 
rate  of  lapinKeaent  was  found  bv  integrating  fquation  (8) 
between  the  following  linits'  Vj,  v^,  and 

Vj, 

i • (9) 


Equation  (^)  was  used  to  nomallre  Tquat  ion  (8), 


d^ 

"n 


- l^r-rV.c  ^ ‘ ^ dV,dV,,dV, 


(10) 


Finally  fqi.ation  flP)  was  intepratcvl  over  all  p')s?»hle 
values  of  Vj  and  V,  yielding  the  pTobahility  density  dis- 
tribution of  the  non  • dimensional  nof'al  vclo^ilv  cof^.ponent  , 
’'i 


Vv^  / 


m ' at 

A plot  of  Fquatior  (11)  is  given  in  Figure  4. 
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Flsurc  i.  Probtbility  Density  Distribution  of  th«  N6n*0iMnitioR«l  NorMl 
Velocity  Coarponent  «C  « ball 


W«1 1 Inffctton  Phtnotena 

The  Bodvl  required  a wall  adher«*d  fluorine  aton  con- 
centration. To  this  end  the  bonding  potential  was 
eaaained.  According  to  Ewiett  (Ref  6)  there  are  two  tjrpes 
of  wall  bonding  phenoaena  divided  into  two  tenperature 
reglnes  in  which  the  atechanisns  are  predoninant.  These 
reglaes  art  known  as  the  low  tenperature  regine,  doninattd 
by  Van  der  Maal  or  physical  bonding,  and  the  high  tenpfra- 
ture  regine  doninaied  hy  chenical  or  activated  bonding 
Because  of  the  operating  tenperature*  of  the  chenical 
lasers  of  concern,  this  study  was  confined  to  the  so 
called  low  tenperature  regine  which  linited  its  applicabil- 
ity to  tenperstures  helow  approkinatelv  BOO  K.  Above  this 
tenperature  the  actual  transition  tc  chenical  bonding  is 
net  kr.;.wr.  for  fi;;crir,r  rr.i  vil!  r.ct  disci.s;cd  l.er^ 


Potential  he \ 1 for  Physical  Adsorption.  A sinple 
nodel  was  cor’*tructed  for  the  potential  well  near  a sur- 
face based  upon  an  adsorption  (attractive)  potential  and 
a repulsive  potential  a<  suggested  bv  l.rnnard • .’ones 
(Ref  IS).  The  nodel  is  described  u'ting  copper  as  an 
exanple.  The  closest  distance  of  approach  of  copper  atoms 
In  nctal  Is  .\!>4  A,  while  that  for  fluorine  atons  is 
l.»’R  A (Ref  16).  The  cquilibriuo  distance,  Rj. , nust  then 
be  the  wean  rf  the  two  distances  or  l.?l  A.  At  this 
equilibriiin  dist.-incc  the  repulsive  forces  are  fron  1/S  to 


1/Z  that  due  to  the  attractive  potential.  Thu»  by  nodellng 
the  repulsive  force  as  40  percent  that  of  the  attractive 
force,  a well  depth  was  arrived  at  by  algebraically 
twBMint  the  two  potentials.  Therefore  an  attractive  potbn< 
tlal.  assuaed  tu  be  V(r).  yielded  a well  depth.  AF^.  at 
the  equillbrlun  position  of 

AE^(hj)  • 0.6  b(Rg)  C12) 

No  reference  Is  aade  to  ar«  activation  energy  and  in  fee* 
no  activation  energy  it  eaperiaental ly  observed  for  physi- 
cal adsorption  (Ref  6).  The  potential  well  is  shown  in 
Figure  S. 

Severe!  authors  have  described  the  functional 
relationship  of  W(r).  The  earliest  of  these  descriptions 
(Ref  IS)  was  the  siwplest  and  yielded  fair  results.  As 
reported  by  several  authors  (Ref  17,]g,  and  IV),  this 
early  node!  served  #•  an  upper  llnit  on  the  physical 
attractive  force  that  one  nisht  expect  to  find.  (>ne  of 
the  wore  recent  papers  dealing  with  the  subiect  was  that 
of  C.  Navroyannis  (Ref  20),  In  his  paper.  Mavroyannis 
derived  a formula  for  N(r)  and  compared  the  results  of  his 
relationship  with  four  previously  derived  relationships, 
those  of  Lennard  - .tones , Rai'dcc)i,  Margenau  and  Pollard, 
Prosen  and  Sachs  (Ref  l!«,n,V8,  and  19).  The  results  of 
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Ph>^ic»l  Adj-irptton  Potential  Well 


his  comparison  showed  that  his  model  compared  as  well  as, 

and  in  most  cases  better  than  the  previous  models  with 

data  for  several  atom-surface  systems.  In  addition,  the 

Mavroyannis  model  made  use  of  readily  available  material 

properties.  For  these  reasons  the  Mavroyannis  potential 

was  adopted  for  use  in  the  fluorine  recombination  model. 

The  formula  for  this  potential  is  given  below 

■hVv 


W(r) 


(13) 


where  r is  the  distance  of  separation  as  in  Equation  (12), 
2 

y is  the  sum  of  the  electronic  charge  of  each  electron 

in  the  adhered  atom  times  the  e.xpectation  value  of  its 

orbital  ra  ius  squared,  “hV.'  is  the  work  function  of  the 

P 

metal  surface  and  N is  the  number  of  electrons  in  the 
adhered  atom. 


The  Interaction  Between  the  Gas  Phase  .-Xtom  and  the 
Wall  in  the  Presence  of  an  .Adsorption  Potential . Given  the 
functional  form  of  the  Mavroyannis  adsorption  potential 
and  the  Lennard- Jones  functional  relationshir-  for  the 
potential  well  depth,  the  well  depth  was  used  to  determine 
the  probability  of  an  atom  of  some  given  energy  adhering  to 
the  wall.  Out  of  several  possible  methods  of  addressing 
the  atom-surface  collision  process  in  the  presence  of  a 
wall  potential-  reported  on  the  literature,  the  soft  cube 
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•odel  pre5(>r)ted  by  Lo|t«n  and  Keck  wa<i>  chosen  (Ref  21). 

S>«e  results  of  this-  node!  were  recently  presented  by 
Fagkii  and  Keck  (Ref  22)  and  Modak  and  Pagni  (Ref  23). 

The  Soft  Cube  Mode  1 . Ax  dictated  by  the  soft  cube 
■odel,  the  surface  of  the  wall  was  pictured  as  atade  up  of 
wall  atons  vibrating  at  a frequency  consistent  with  the 
bulk  Pfbye  tenperature  in  a direction  nornal  to  the  plane 
of  the  surface.  The  overall  picture  was  that  of  an 
Infinite  checkerboard  with  each  square  made  up  of  one  aton 
vibrating  randonly  into  and  out  of  the  surface  at  the 
[iebye  frequency,  w^.  During  a surface  interaction  only  the 
velocity  conponent  of  the  gas  phase  ato*  nornal  to  the 
surface  was  considered.  The  use  of  the  noraal  velooty 
coMportent  ruableu  data  correlation  which  was  not  able  to 
he  aiade  using  the  total  velocity  in  the  case  of  low  inci- 
dent angle  a t om  • sx»r  f .see  collisions.  The  gas  phase  atom 
was  furthc!  assumed  to  interact  vith  nnlN'  one  surface  *to*. 
The  variovis  frequencies  and  potentials  were  then  modeled 
as  shown  in  1 igure  (• . The  sh.ipr  of  the  well  hevond  the 
interaction  posit  lou  was  of  n«*  importanc  e,  ^he  gas  atom 
was  simp))'  given  additional  velocity  above  its  Bolt.rmann 
predicted  g.s«  phase  vrlocitii  consistent  with  the  well 
depth.  As  dlscusset!  bv  'fodai  ami  Tagni  the  interaction 
spring  constant,  was  mridelrd  by  an  exponential 
repulsive  potential,  V(r), 
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Figure  6,  Collision  Model 


V(r)  • C I I ’'e  • 


«h«r«  r is  the  position,  « is  the  distance  over  which  the 

interaction  took  place,  and  C and  t ar*  constants  adjusted 

to  naten  experinental  data.  The  spring  constant  was 

obtained  from  fquation  (14)  by  dl f fetent iat ion  to  o tain 

the  force,  and  expansion  of  this  force  around  r • D keep 

inf  on!;,  up  to  the  linear  tern,  Matching,  the  boundary 

condition  that  the  force  " 0 at  r • 0 yielded  a spring 

constant,  k , 
f 


k * 

K 


.2  !>' 


where  1>  is  the  well  depth . The  '.pTing  cc*nstant  for  the  sur 
face  .ston  w;i\  givnt  hv 


ri 

S s 


A colli'iori  heg.»n  vhrn  .i  ph.»‘-p  , v ire  in 

contact  *«  i t h the  •'oviuc  end  of  thr  vprnic  V^,  ‘ ron  ithic 

point  until  the  volli'-lon  onded  the  v.  trr  at  led  uv  u 
sioplc  und.mj'Cd  tvo  oii".-,  •iv'.trn.  A volli'^Joi  then  tool 
plavc  vit>i  a ,;nN  5'h.sse  ati>r  of  velocity  v 


(H) 


where  is  the  Foltzaann  predicted  random  surface  normal 
component  of  velocity  (see  Equation  (11)).  The  Y(0)  posi- 
tion of  the  spring  at  collision  was  also  random  but  equal 
to  2(0)  The  initial  conditions  were 


2(0)  - y(o) 

(l«) 

V > V(0) 

(19) 

Equation  (I'l)  is  interpreted  to  tncar.  that  if  the  velocity 

t 

V were  less  than  V(C)  there  would  have  been  no  collision. 

The  collision  ended  at  when  Y(t^)  ■ l(t^, ),  At 

% 

that  point,  if  1/:  m (v(t  ))*'  F (R,-)  the  atom  vas  fapped 

1^  c *•  j 

and  if  1/2  a.  (v(t  1,  the  alow  escaped  the  well 

f c ■ w I, 

and  avoided  trap}<ui};. 

Referring  to  I i pure  f' , the  equations  ot  wot  ton  for  the 
collision  j’iDte*is  ave  given  by 


g s g 

mV-k2*kV"0 

g ? g 


(20) 
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Tit*  detailed  solution  of  Equaitions  (2H)  is  given  in 
Appendix  A.  The  results  of  this  solution  yielded  the  func- 
tions t(t),  2(t)  and  Y(t).  Using  the  collision  criteria 
above,  stick  or  no  stick  decisions  were  atade  for  given 
sets  of  initial  conditions. 

Kate  Process  am?.  Equations 

The  stick  or  no  stick  decisions,  as  described  in 
Appendix  A,  required  40  calculations  for  each  of  26  initial 
non-dinensional  velocity  conditions  predicted  by  a 
Boltzmann  distributed  gas  for  a given  temperature.  Twenty 
temperatures  were  run  between  0 K and  flOO  K for  each  gas- 
■etal  system.  The  result  of  this  process  yield«'d  an 
initial  or  clean  surface  sticking  coefficient,  S^(T), 
defined  as  the  statistical  orobabilitv  that  any  given 
collision  ol  a gas  phase  atom  with  the  surface  at  tempera- 
ture T would  result  in  the  atom  h**ing  trapped  by  the  sur- 
face. The  actual  scickinc  coefficient  as  a function  of 
surface  coverage,  ((»  ■ I being  a fo.My  covered  ronolayej 
of  adhered  ruri.icc  atonv),  has  exper ir.ent a 11  y verified  by 
Christmann  et.al,  (hef  241  to  be 

S(T.P)  - (l-f-)S^n)  (21) 

The  sticking  coefficient  was  then  used  in  describing 
the  rate  equation  for  the  time  rate  of  change  of  the 
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surface  adhered  atom  concentration  bv 


oAp*  «p*  - 


(22) 


where  n.  is  the  surface  concentration  and  K.  is  the  rate  of 
inplngenent  of  F atoms  as  described  earlier.  The  depletion 
of  the  surface  concentration  described  by  the  last  two 
terms  of  Equation  (72)  was  given  by  the  thermal  desorption 
rate,  tines  a surface  coverage  9 (hef  6)  and  the  detach* 
went  caused  by  reconbinat ions  was  represented  by  a sterlc 
factor  S^.  the  surface  coverage  0,  and  themal 

desorption  rate  used  was  that  given  by  Classtone.  laidler 
and  Eyring  (Ref  25) 


6 


kT 

'•a'Ti 


FT 


(25) 


where  was  the  number  of  surface  sites  per  unit  area 
(netal  surfase  atow  packing),  h the  Planck  constant,  and  F.|^ 
a desouption  energy.  The  packing  used  was  that  of  the  «ace 
centered  cubic  for  all  surfaces  considcjcd  (Ref  26).  The 
steric  factor  is  the  ratio  of  the  gas  phase- surface 
adhered  atom  collisions  resulting  in  a recoaibinat  ion  to  the 
actual  number  of  these  collisions. 

The  surface  coverage  of  interest  in  Fquat ion  (22)  was 
of  course  that  of  the  fluorine  atoms.  Since  the  actual 
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nuaber  adhered  atoms  depended  on  the  nuabrr  of  available 
surface  sites,  the  presence  of  adhered  atoas  other  than 
fluorine  effectively  decreased  the  surface  area  available 
for  fluorine  to  adhere.  In  this  light.  Equation  (22)  was 
rewritten  as  follows 

d(n^L 

— ' SjApbl.  (24) 

where  9^  represented  the  fluorine  surface  coverage  and  9^ 
represented  the  total  surface  coverage  which  included  ether 
as  well  as  fluorine  atoas. 

To  quantify  9^  a rate  equation  was  written  to  include 
the  effect  of  any  diluent  which  night  have  been  present. 

This  rate  equation  was  given  by 


where  f-quation  reflects  only  one  subtracting  term, 

that  of  thnnAl  desorption,  for  one  diluent  then, 

9j  “ ••  Op.  assuning  molecular  species  do  not  adhere. 

Solving  Equation  (24)  and  (2S)  for  the  steady  state 
condition  yielded 


(SJ 


(1 


(29) 
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•f 

* * t?;Pf  *F  * tv; 


fl 


•o> 


(H) 


Further,  tolvinf  Equations  (26)  and  (27)  tinultaneously 


‘ ■ <*D>o 

"■^ww 


(Zl) 


uhere 


<«F>o  = n sr~  <”) 

' • ts;7;  V ts;t; 

f'n'o  • -t; <»») 

(Op)^  and  vcrc  defined  to  correspond  to  the  actual 

coverages  of  one  in  the  absence  of  the  other. 

Having  soiveu  for  the  surface  coverage  of  atonic 
fluorine,  the  reconbinat  iui«  rate,  Vj^,  was  given  sinply  ttf 


V 


R 


(31) 


Finally  the  reconh i nat ion  coefficient,  v,  defined  as  the 
ratio  of  fluorine  ston  wall  collisions  icsulting  in  a 
reconhinat ion  to  the  total  number  of  fluorine  aton  wall 
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col  Hu- tons  , was  }jivc'ri  by 


* ^ ■ ^^0**  Presnure  Pipe  Flow  with  Wall 

Kecoabinet ion 

A sinple  Model  of  one-di»enslon«)  voluvetric  fluorine 
depletion  applicable  In  low  pressure  pipe  flow  w(  s 
developed  to  siaiulate  the  Ultec  cxperlaient  (Ref  15). 

Ultee  ExperiMent 

The  Ultee  data  were  taken  in  the  experivental  apparatus 
described  pictor ' al ly  in  Fleure  ?.  Gas  of  known  F2  concen- 
tration and  a diluent  far|>on)  flowt-d  froai  thr  fas  source 
into  a Microwave  di^charfc  cavity  where  the  Molecular  fluo- 
rine was  dissociated  to  atonlv  fluorine.  Mcasurenents 
taken  at  the  exit  of  the  cavity  without  the  test  sauple 
tube  attached  showed  that  approxiaatcl)  lb  pe/cent  of  the 
Fj  was  dissociated.  This  fixed  the  jas  Mixture  of  F,  F^ 
and  Ai  at  the  cntr.ince  of  a “5  cm  X 1.1  cm  (I.  D.)  test 
saMple  tube  wade  of  various  catalytic  nctals.  interest 
here  was  the  nickel  sunple.  The  pressure  drop  across  the 
saMpIr,  thr  eqiiilihriu*  trnpcratuic  of  the  sa»p\e  and  the 
flow  rate  of  the  ras  were  isonitMrcd  and  recorded.  An  I SR 
spectroMcter  located  at  the  end  of  the  test  sanplr  yielded 
inforwation  fro»  which  F ami  Fj  concentrations  were  calcu- 
lated. The  fas  then  left  the  test  apparatus  through  the 
puMping  svste»,  For  each  data  point  the  syster  was 
allowed  to  come  to  equi 1 ibr inn  at  the  particular 
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Fiiture  Ulte<?  Expcriatntcl  App«r«tus 


t««pcrature  of  interest,  always  Ra^pini  a conMant  initial 
praasurt  and  flow  rata  throufhout  tha  axpeTlaanlrJ  range  of 
taaparatures . 

Analytic  Mode? 

Undar  the  assuaiption  of  a one -dimensional  dapendanca  of 
tha  concentration  on  axial  position  in  the  test  saaipla,  the 
equation  of  the  process  was  written  as  follows  (see  Figure 

• ) 


dn  _ 4yn  / kT 

Jx  " Or  V Tim  j 


(531 


where  n is  the  nuaiecr  densitv,,  * tite  axial  position,  r the 
tube  radius,  t the  recombination  coefficient,  U the  flow 
velocity,  and  n (kT/^i'n)  the  rat?  of  impingement. 

F.quatio;i  (3^)  has  an  analytic  solution  under  the  assumptions 
of  a constant  rc«.onbinat  ion  coefficient  r and  an  initial 
number  densitv  n 


4>* 


(54) 


Fquation  (34)  can  be  uved  to  dctc laine  an  overall  effective 
recombination  ccefficlent. 

Since  the  recomhin»t ion  coefficient,  as  described  in 
Section  IT,  depended  on  various  parameters  vhich  changed 
with  position  along  the  test  sample,  more  accurate 
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igure  I,  One-Diaenilor.®!  Pip«  Flow  Conc«ntr«t ion  Doplotion 


simulation  of  the  actual  experiment  was  made  by  a step  by 
step  solution  of  Equation  (34)  for  an  axial  grid  spacing  of 
Ax. 


n(x^  + Ax) 

HTxTl 


e 


4y(x-)Ax 

FU 


(35) 


where  the  values  of  y(x)  and  n(x)  were  calculated  at  each 


station  to  predict  the  change  in  the  number  density  to  the 
next  station. 


Applicability  Criteria 

As  a measure  of  the  applicability  of  the  sole  assumed 
dependence  of  the  fluorine  concentration  on  axial  position, 
the  characteristic  times  for  diffusion  and  depletion  can  be 
compared.  A typical  time  for  diffusion  to  negate  gradients 
can  be  represented  by  the  relaxation  time  for  diffusive 
equilibrium  in  an  enclosure  of  dimensions  r (Ref  27) 


^DIFF 


(36) 


where  V is  the  diffusivity,  X is  the  mean  free  path  and  v 
is  the  average  molecular  velocity.  This  time  can  be  com- 
pared to  the  time  necessary  for  destruction  of  all  the  atoms 
in  the  volume  V by  recombination  at  the  wall  neglecting  the 
effects  of  diffusion.  The  depletion  time  is  given  by 
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where  S is  the  wetted  surface,  for  validity  of  the  one 
diaiens ional  assuaptlon  of  the  concentration 

*IIEC  *DIFF 


(3«) 


IV. 


owpaT'i  stjn  of 
with  lKj>(i*rii»ental  Data 

St ickin£  Coefficients 

The  solution  of  the  equations  of  Motion  fo?  the  colli" 
tion  process  were  integrated  over  the  Bolttaann  distribution 
given  in  Figure  4 to  obtain  the  sticking  coefficient,  S^, 
for  fluorine,  argon,  and  heliu*  for  three  surfaces  represen- 
tative of  fluorine  laser  no::le  Materials;  gold,  nickel  and 
copper.  The  results  are  given  in  Figures  9,  10.  and  11. 
CoMparison  of  the  figures  in  light  of  the  effect  pf  diluents 
Of*  re  ronhinat  ion  rates  (Fquat  ions  ('81  through  (Tl))  *s 
discussed  earlier  indicatcu  that  argon,  when  used  as  a 
diluent,  should  Inivc  an  appreciahlt  effect  on  the  rcconhin-i' 
tion  rale  of  fluorine,  while  the  hrl'.uai  prohahl>  would  not. 
Furtho  Bcniinr.  of  the  hellus-  ihluciit  svstre  is  wade  later. 

tUtce  (ReC  Itl  rc|H*r*  ’t!  i . a slightly  higher  reson 
blnatior.  rate  J-,avo  octuired  or  coiipcr  than  I'n  nui^el. 

I xa^’in.it  ici!  of  fi.'uc  > fun. th.it  tfiis  i :«  i on" :-tcnt  with 
a slightlv  higfri  ••tivriUf.  v "i*  f f , v i rut  r«(  t):i  fluriuu'  ur» 
copper. 

(,*o»y>a  i i "-or,  y>f  Fij  c Fly’‘  C.i  K u 1 at  j on--  with  I'.vt  • 

The  pftVMv.t)  i‘aiarct»r‘i  of  tfu  UlUe  isjM  iucut  .«ir 
given  iti  Table  I .and  the  ( oeput  .it  ) rr^.t  1 p a r a**  f ’ r i :*rr  given 
in  latfSe  II.  IKing  the  app  I n,  .‘b  i 1 1 1 v i»Ueii-"i«  fiur  ! qua 
tion  (^8)  It  was  found  th.it  ca  Kul.it  i<*n<  b.ised  on  tne 


;?nt  for  Atonic  Fluorine 


Sttcking  Cpcfficient  foT  ArfOn 


one'dinen5lonal  assumption  for  tho  conditions  of  the 

Ultee  experiment  were  applicable  for  a recombination  coef- 
ficient  as  high  as  10  Thus  the  one*diinenslonal  assump- 
tion vas  good  for  all  available  data  points  as  will  be 
aho«m. 


Table  I.  Physical  Parameters  of  Ultce  Experiment 


PARA.HF.il  R 

VALUE 

PRESSURi. 

Z . 7 TORR 

rtoif  vuonTv 

180  cti/sec 

MOLE  FRLiVCTSON  A7 
SAHPU  ISTRAVn 

(H 

l»2,l 

|AtJ 

0.01 

0.0^ 

0.90 

TEMPI  RATI, 'H!  KASt  l 

0 *(  to  400  *C 

SAMPU  WnRlAl 

hlCFEL 

SAMPM  i,r<„ 

1.1  cm 

romput at  3 jjn,'*  based  on  T^ivat  lOJl^  ami  (IS)  were  per’ 

forned  for  cn-’ i or;  wjitb  the  expo  3"  ; a I data. 

The  I e Milts  of  !l'-ese  cor  p>jt  at  3 on^  wre  pu‘ 'tented  'im  5'3};me 
12.  Shown  in  Heme  Id  is  the  "rri'cnt  f reffijun jnp  at  t?ie 
end  of  the  niiKei  sar;jile  tube  over  .i  rar)j.e  of  ^tead>  state 
tcEju ')  aluros  . The  only  paraneleis  viijth  were  arbitral  ily 
cho'.rt!  in  ordm  to  have  fJic  rnirdol  natch  the  data  were  the 
deserptiew  ''rxrp.es  ^*'I?\\r 


, and  the  factor 


(sec  Table  II)  , The  desorption  rnergltris  required  to 
match  the  data  were  found  to  he  siailar  in  ae^nifude  to 
those  for  other  gas -surface  s/steas  reported  by  Eaaett 
(Ref  b).  Further,  the  choice  of  the  desorption  energies 
uas  rather  clear  since  a higher  or  lower  energy  caused  the 
curve  intersection  with  the  tenperature  axis  to  shift  to 
the  right  or  left,  respectively.  Also  as  an  additional 
Indicator  of  the  correct  choice  of  desorption  energy,  a 
higher  or  lower  energy  caused  a shallower  or  ste%?er  slope, 
respectively,  although  this  effect  was  not  as  pronounced 
as  the  intersection  shifting.  The  effect  of  varying  the 
steric  factor,  S^,  is  shown  in  Figuje  15.  A steric  factor 
of  0.1  seened  to  best  fit  the  data. 


Table  II. 
rAR-AWTIR 

Ar 


' * Ii^Ar 


f^^Ar 


Conputat ional  rarameters 
VALUI 

6.632S8  X IP  g 

3.1S451  X 10*^'  g I 

I 

.5.587S  X 10“' 

.j.  I 

1.^4:  X 10  '•  erg  j 
1.175  X 10'* ’ erg  I 
P.  1 

Flgr,.e  9. 

Figure  10. 
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( ''c.  ) 

Hfmxc  n.  Ifffci  of  Vftiyin^:  Stfru  lactoi 


<s* 


By  using  tfie  nodvi , regions  of  the  date  were  Inter- 
preted and  labeled  as  shown  in  Figure  12.  Region  A (0  *C 
to  40  *C)  is  that  region  in  which  the  'rail  reaction  was 
dominated  by  the  presence  of  adhered  argon  atons 
'‘poisoning"  the  reconbination  process  by  not  allowing  free 
surface  sites  for  the  fluorine  to  attach,  by  use  of  liqua- 
tion (30)  It  can  be  seen  that  prior  to  reaching  soaie  acti- 
vation tenpersture  for  thernal  desorption  (l.e., 

4p  * < argon  copspletely  covered  the  suiface 

therefore  rendering  cqu  .l  to  :,ero.  As  the  tenperature 
rose  the  therwal  desorption  becane  inportant  and  the  argon 
began  to  allow  portions  of  the  surface  to  be  left  vacant  so 

the  f Uie  could  attach  and  subsequently  reconbine. 

In  region  A,  (40  ♦<*  to  195  *0  w,is  ins i gni  f U' ant  as 
coapafed  to  thus  the  surface  coverage  was  deter- 

wined  prinarily  by  5^  and  so  that  rquaiion  (29)  could 

be  approx  mated  by 

• ‘y-—  (J9) 

‘ * rs„T, 

Finally  m region  r (195  *C  to  -tOC  *n  the  desorption 
of  fluoriiir  was  activated  and  e\rntu.sny  doninaicd  the  wall 
re.y<  t ion . 

At  this  point  it  is  of  interest  to  aention  Figure  14. 
r.quation  (54)  was  vised  in  coniuruti  n with  the  UStec  <".Tita 
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to  obtain  an  avera|(e  Mall  recoobinet ion  coefficient.  The 
result  of  this  coMputation  is  given  in  Figure  14.  The 
unknoini  regi'nn  rtfars  to  the  teaiperatures  at  which  no  F was 
detected  at  the  exit  of  the  test  sasiplc.  Thus,  at  an  aeer- 
age  Y of  approxiaately  10  all  the  detectable  F had  been 
reoKsved  from  the  flow  by  wall  reconbinat ion.  This  seeas  to 
indicate  that  within  the  range  of  the  expenaental  data  the 
assuaptlon  of  one-dincnsional  depletion  was  applicable. 
Equations  (26)  through  i30).  however,  show  the  dependence 
of  Y on  the  local  nuaber  density  in  the  terns.  The 
aaxiaua  recoabinat ion  took  place  at  the  entrance  of  the 
test  saaple  where  was  a aauiaua.  Figure  IS  shows  the 
teaperatiire  dependence  of  y based  on  the  naxiaun  nuaber 
density.  Using  the  aaxiaua  y criteria  the  experiaental 
data  range  still  falls  within  the  one  diaensional  assuaptlon. 
More  iaportant,  however,  is  the  capability  of  dctcraining 
the  actual  teaperature  and  nuaber  densiiv  dependence  of  y 
using  the  aodcl.  Without  the  nodel  the  best  one  could  do 
is  to  dctc"i®ine  an  average  y from  Figure  14  extrapolated  to 
a aaxinuR.  Clear 1>  the  actual  «axipua  is  nearly  Two  ordcis 
of  Kagnitude  higher  than  would  he  estinnted  w.thout  The  use 
of  the  node].  This  aistaken  naxir.ur,  i^  in  fact  vhal  was 
reported  by  Ulter  (Ref  l.\).  A set  of  naxinun  reverb i nat  ion 
coef  fit- lent  curves  for  several  gas  pressures  for  90  percent 
dissociated  <"luorine  is  given  in  Appendix  While  the 


4S 


X 

CMXUATcD  f£.SUL' 

\ 

\ 

TOR  mm 

% 

% 

viujon 

\ 

S 

CALCUUltD  RESU. 

R»  ICLIIM 

/ V 

mUJENT 

/ \ - 

C^CLLATED  F«0H 

IXTIWP0LAT1ED 

ivixrui 

ICFEfCirf  j; 
ORTA  («W  U ) 

TFJf  iKAlU*  i K ) 


F»fur<‘  IS.  V.’ixc'su’*  riuotirv  KjH  f'.i,  x itJ:  init  k'ii  ( cif  f »i.  nnt 
vs.  ■!  rspr  r.u  mj  c for  a Prvs'^jrc  of  2."  lor 

901  Arfon  ial.Kn?  Jlovi,  .iiiJ  101  fl»Hn"»nc 
Dissov  rat  ion 


curves  for  90  percent  dissociation  are  different  than  vould 
be  predicted  by  10  percent  dissociation  the  characteristic 
trend  for  increased  total  and  partial  pressures  it  the 


saae . 

based  upon  Figure  H,  it  is  prcdUteJ  that  in  a 
siallar  test  run  with  a heliuai  diluent  region  A would 
either  be  nonexistent  or  greatly  reduced  fro«  that  in  Fig- 
ure 12,  while  regions  B and  C would  show  no  change.  The 
•axinun  rrcoeh^nat ion  coefficient  with  a heliun  diluent 
should  roughly  correspond  to  the  curve  in  Figure  IS 
represented  by  the  dashed  line. 
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Effect*  of  Recowbir.at  iion  in 
Cheaicel  Laser  Xo ; ; 1 e s 

Flow  Equat ions 

A two  dlnensionel.  steady  state,  laoinar,  explicit 
boundary  layer  rode  was  con  tructed  to  approximate  nostle 
flow.  The  code  was  constructed  for  the  purpose  of  denon- 
strating  the  applicability  of  the  w-!l  recombination  model 
to  real  flow  system*.  The  governing  equation*  are  as 
follows  under  the  usual  boundary  layer  assumption*: 
COKTIKUITY: 


MOMfXTifV 


SPf  CM  S 


jy(ou)  • ® 


(40) 


pujj^  ‘•'SV-  Jx 


(41) 


»K 

puji- 


if 


■ . h * U l'i) 

) iyV)iv  / 


(45) 


' Ue**" 

p( 


* 


Jp 

'Mt 


t*’  i. 


91 


*v 


^ W K * ^ V , «■■  ♦ 


T ‘ 


t ' ? I . 

i , . 


(43) 


4S 


wh«re  X Is  the  axial  position,  y the  transverse  position 
(see  Figure  16),  r the  density,  ti  the  x ■ conponent  of  veloc- 
ity, V the  y-roaponent  of  velocity,  P the  pressure,  u the 
viscosity  of  the  Mixtuic,  the  rate  of  production  of 
species  j,  the  diffusivity  of  specie  j into  the  Mixture, 
the  Bole-aass  ratio  of  specie  j,  the  specific  heat  at 
constant  pressure  of  the  Mixture,  h^^  the  enthalpy  of  specie 
I,  and  fc  the  theraal  conductivity  of  the  ■ixture.  Von  Mises 
trans fornat  ion  for  coapres^iblr  flov  was  used  to  s.iaplify 
the  equations  for  insertion  into  the  code.  Introducing  the 
streaa  function 
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the  difference  equaticn*^  /or  I'tju.it  ions  (4C),  (4(«)  , and  (4’’) 
are  iricluded  in  Appendix  ' . 
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Thermodynamic  Properties 

The  enthalpy  and  specific  heat  at  constant  pressure 
for  the  individual  species  were  calculated  in  the  form  of  a 
quintic  polynomial  with  temperature  as  the  independent 
variable  as  follows: 
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where  the  coefficients  a,  b,  c,  d,  e,  and  f were  taken  from 
Reference  2,  and  R is  the  universal  gas  constant.  The 
density  was  calculated  by 
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Homogeneous  Reactions 

Two  free  stream  reactions  were  considered,  that  of  the 
dissociation  of  fluorine  and  that  of  the  recombination  of 
fluorine . 
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RECOMBINATION 
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where  M represents  any  second  or  third  body,  respectively. 
Here  the  dissociation  was  defined  as  the  forward  reaction. 
The  forward  reaction  rate  used  was 


Kj  = 4.2  10-11 


(53) 


taken  from  Reference  5.  The  backward  reaction  rate  used  was 
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taken  from  Reference  28. 
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Transport  Properties 

The  transport  properties  for  the  flow  constituents  were 
taken  from  Reference  29. 
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where  Oy  is  the  reduced  collision  integral  curve  fit  as  a 

function  of  kT/e^  and  is  the  molecular  weight  of  species 

i.  The  values  for  o.,  the  cross  section,  and  c.,  the 
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maximum  energy  of  attraction  between  colliding  molecules, 
were  taken  from  Reference  30.  The  viscosity  for  the  mix- 
ture was  calculated  using 
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where  is  the  mole  fraction  of  species  i 
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THERMAL  CONDUCTIVITY 

The  thermal  conductivity  of  the  flow  constituents  was 
calculated  using  the  Euckcn  approximation 
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The  thermal  conductivity  of  the  mixture  was  calculated  by 
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is  again  given  by  Equation  (57) . 
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DlFFUSrVITY 

The  diffusivity  was  calculated  by  first  calculating 
the  binary  diffusion  coefficient  by 
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ttp  was  curve  fit  as  a function  of  kT/c^^  where 
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and  O-j  was  given  h> 
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with  siiiiiar  equittions  for  and  Pj^  (Ref  51) 

More  Than  Three  Species 

For  BK>re  than  5 species  only  trace  species  were 
handled.  KWcn  wore  tnan  3 species  were  ptesrot. 
the  first  5 were  (iven  by  Equation  (64)  and  the 
reiaaining  trice  species  by 
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Roundary  i t ions 

The  no  slip  wall  condition  for  velocity  was  used.  The 
tcaperiturc  at  the  wall  was  specified  either  as  sone  con- 
stant value  or  as  a ^u»ctior^  of  position  alonjj;  the  wall. 

The  wall  con/.i  t ion?  for  v.  once  nl  rat  i on  of  atonic  and  molec- 
ular fluorine  were  giver.  <is  fellows  a function  of  the 
recoTObinat  ) PI.  lorfficicnt  (deiivation  included  in  Appendix 
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The  domain  of  the  calculation  was  from  an  axial  posi- 
tion* ^ ® 0.0  at  the  nozzle  entrance  to  x = 0.4  at  the 
nozzle  exit  (see  Figure  16) . 

Results  of  Nozzle  Flow  Computations 

The  boundary  layer  code  was  run  for  predictions  of  the 
effects  of  recombination  on  laser  nozzle  flow.  The  pres- 
sure profile  and  nozzle  contour  for  a representative  slit 
nozzle  are  given  in  Figure  18  (taken  from  Ref  4) . While 
not  the  most  commonly  used  slit  nozzles,  it  was  similar  in 
size  and  basic  contour.  The  input  conditions  for  these 
computations  are  given  in  Table  III. 

Computations  were  run  for  three  wall  catalytic  condi- 
tions: (1)  noncatalytic  wall,  y = 0.0,  (2)  fully  cata- 

lytic, Y = 1.0  (Ferrell  et.al.  conditions)  and  (3)  step  by 
step  computation  of  y based  upon  the  recombination  model. 
The  results  are  presented  in  Figures  19  through  28.  These 
results,  as  mentioned  previously,  were  based  on  the  assump- 
tion of  laminar  flov/.  The  highest  Reynolds  number  based  on 
the  momentum  thickness  was  approximately  65.  This  was  well 
within  the  laminar  regime  as  transition  to  turbulence 
occurs  at  around  360  (Ref  32) . 

Velocity.  Figures  19  and  20  show  the  effect  the 
recombination  had  on  the  velocity  profile  at  the  throat  and 
exit  plane  respectively.  As  can  be  seen  there  was  no 


Figure  19.  Velocity  Profile  at  Throat 


major  difference  in  the  velocity  profiles,  only  a slight 
lessening  of  the  velocity  gradient  at  the  wall  with  recom- 
bination, decreasing  the  most  for  the  case  where  y = 1.0. 

Temperature . Figures  21  and  22  show  the  effect  the 
recombiantion  had  on  the  temperature  profiles  at  the  throat 
and  exit  plane,  respectively.  The  throat  profile  showed 
very  little  change  with  only  a slight  steepening  in  the 
temperature  gradient  at  the  wall  with  increasing  recombina- 
tion. The  exit  plane  profile,  however,  showed  a more 
dramatic  effect  both  in  an  increasing  slope  at  the  wall  and 
a shifting  of  the  profile  near  the  boundary  layer  mid- 
portion . 

The  overall  wall  effect  appeared  to  be  a slightly 
increased  heat  transfer  to  the  wall  due  to  normal  conduc- 
tive heat  transfer  (see  Figure  22)  . To  this  sliould  be 
added  the  heat  released  in  the  recombination  reaction.  For 
this  case  the  temperature  was  assumed  constant  along  the 
wall  as  indicated  in  Table  III,  but  an  actual  determination 
of  the  nozzle  wall  temperature  profile  should  not  neglect 
the  role  of  the  energy  of  the  reaction.  This  energy  is 
dumped  directly  into  the  v/all  as  indicated  by  Meyerson 
(Ref  33)  prior  to  being  either  conductively  carried  away  by 
the  wall  or  recoupled  into  the  flow  througli  normal  heat 
transfer  processes.  This  additional  heat  transfer  turns 
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out  to  significant,  vjp  to  thrc«  tines  the  nornal  con- 
duction as  can  be  seen  in  Fijiurc  ’5. 


Table  III.  Input  Conditions 
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cavity  co:  j'utat  . The  typic.ij:l  visctMJS  heating;  dis- 
tort ten  of  the  pvi.1  f 1,  Ic  as  i nd  1 c ;6 1 oiH  I- * the  norc  .ita  lyt  sc 
Ciirvi'  15;  ?'t,-:5.TC  .’J,  js  affcvtcd  hv  the  i!;  f fs;;.  icn  of  5 aj;d 
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hiKhcr  cii t ha  1 j.s-  of  t^o  f .ttOT-i  jis  ctVi'-pa red  to  one  !\  r,o.le- 
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culc  tetidod  to  <’r,:o  esu'rps  i»i:t  of  the  oa  ils  t roa-  thereby 
SMOOThsi'.K  ,niit  the  profile. 
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Ch«'«^i_c_al  i <■  The  effect  of  recombination  on 

the  concent  rat  5 cn  piTfi'cs  was  of  west  interest  and  xs 
exaninecl  in  Fixtures  24  t:irough  2S.  Figures  24  and  TS 
five  the  profiles  at  the  throat  and  exit  plane,  respec- 
tively; figiire  2^  'hows  the  depletion  of  f and  jrowth  of 
Fj  with  axial  post  ion  along  two  strean  functions  in  t'" 
norsJe;  figure  r”  shows  the  loss  of  f,  and  figure  28 
tains  the  reco'''hi  nat  < on  crefficient  as  a function  of  -mial 
position.  X3thcM,!j:f  Figure  “ 24  and  2S  show  a differeiHC 
between  tVt  fulJv  cat*ly!;v  versu*  the  r.odel  , it  is 
notiteaH>  of  s Jiffermce  thaa>  the  one  and  a half 

order  of  raj..rii  t-ide  difference  in  the  rccmtfinat  ion  creffi- 

cicnt*-  i'‘'Cr  rtg;ae  ‘S  . '.’he  icasc;.  for  this  can  he  fouviJ 
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Finally,  rearranging  the  terais  slightly 
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Thus,  the  actual  rate  at  which  F is  depleted  froB  the 
systva  ?s  dependent  on  both  the  diflusivity  and  the  recoa- 
bination  coefficient.  Fxanination  of  the  two  terns  in  the 
denoninator  of  Fquation  ('*0)  yields  the  figures  in  Table 
IV.  The  table  shows  that  while  the  y’s  differed  by  1",  the 
actual  dcpleti>.n  was  liattcd  to  a difference  of  only  S at 
the  exit  and  4 at  the  throat  due  to  only  the  terns  in  the 
denoain.itor  of  Iquaticn  ("0).  Coupled  to  thi<  Itaiting 
influence,  the  actual  F atoir  concent r.it  ton  in  CFU.  2 linited 
still  further  the  difference  in  the  twe  solutions  so  that 
the  final  differences  were  reduced  f r<n  1”  vo  approxinatel y 
i . 4 . This  ratio  is  verified  hy  Figure  2“,  where  the  total 
anoiint  of  atorsc  fJuori.ie  removed  by  rcconhinat  ion  was 
52  percent  c>>t’ vu rr<*  with  4S  ;<rccnl  nneved  when  - • l.C. 


VI  . Conclusions  and  Rccoftisondat  ions 


Tho  present  work  dealt  with  the  construction  of  a 
theoretical  Rkodel  for  the  hcieroiteneous  reconbinat ion 
reaction  of  fluorine  at  a netal  surface.  Usinit.  the  Rtdeal 
■techanisR  the  reaction  was  divided  into  three  distinct 
processes  (1)  the  bundtnc;  process  of  atonic  fluorine  to 
the  wall,  (2)  the  collision  and  reconbinat ion  process  of  a 
pas -phase  fluorine  a tor  and  a bonded  fluorine  aton  and  (M 
the  rate  process  at  vhich  the  honJtuK  and  ret. onb mat  ion 
takes  place  cer,' idet  inf  the  surface  coveraj^e.  Having 
treated  each  of  the  procf*:.*e*  in  sere  detail,  iht'^e  pro- 
ces'sev  vrre  co'^hmed  J.  fi-rr  the  fii  .il  wall  rct'o*hm.it  non 
node  1 . The  eoi’i  i vj«  avajtcJ  to  vri  • J ; - ei)'  . c.._  1 pije 
flow  and  n.  c.-’ini?  ton-  radc  for  vn  vit*'  data  taker, 

by  hr.  fft'j'er  *.  ’.'Iter  j Vinited  Vit^ruft  J^'i  carvh  l.rlcra 
torsos.  '■  mil.  3 1 '»■  tie  w.i‘  j'laco.l  jr,.to  j t wo  • d i P\-r.>  ; c:iu  I, 

stcadv  -'.ate,  '...•iMT  V.  uT...'.>rv  S.-.vnj  ..rdc  t*  d»'T  j i the 

.'lisc  fro~  "rCi.'"  f jj-- t vjijS  rrv-'''T;- 

b,ir!„i!  ; ( tn  • < ■ . Ic  of  t ' v *.  ;i  i r,  f 1 . r ;■  f » n>.  - . . ,.l 

la*.e  . 
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(1i)  a low  icuprtaturc  regiite  (less  5S0  %)  in  which  the 
recoi»bi:iat  ion  reaction  was  poisoned  by  argon  doninat  in^;  the 
furface,  (2)  a nid- temperature  regime  (?50  K to  47S  K)  in 
which  the  argon  was  thermally  desorbed  and  the  fluorine 
thermal  desorption  was  not  significant  enough  to  suhstau* 
tially  decrease  the  surface  coserage  of  fluorine  and  (5)  a 
high  temperature  legime  (fS  K to  .''S  K)  in  which  thermal 
desorption  playvd  a s ignt ficant  role  in  decreasing  the 
fluorine  surface  coverage.  Ifhiie  the  effect  of  other  dilu' 
ents  on  the  ? regme  ^ specuJat  rvc , it  is  predicted, 
based  on  conputat  cuntairci  htmn,  that  the  Ttco.nbina- 
non  cvrffu'scM  would  * t'  »•  ifni  f uMnt  li<  hightt  in  the  low 
te*>pt  r at  r rcf.'"  .y  h<lrje.  dj  Suont  , The  third  regiar, 
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lun  by  Ultt'c  should!  ho  ruri  hiitK  a hcliun  diluent  to  fully 

vorxfy  thr  insignificance  of  attarhatent.  Additionally, 
future  efforts  should  be  addressed  to  the  prohles  of 
interfacing  this  nodcl  4rlth  production  flow  codes  desifned 
specifically  for  analysis  of  ch^aical  laser  nozsles  such  as 
the  one  described  in  Reference  2. 


Hblto£r>f.hy 


12. 


U. 


14. 


IS. 


1ft. 


17. 


IS. 


19. 


20. 


2 3 .. 

22. 


25. 


Valoinrr.,  ft'.,  Sirinp,  B.,  anJ  Maclean,  P.  I.  "Veasure- 
nent  of  Flu''>rinc  \tcif  Concent  rai  iuns  and  «leco»bina- 
tlon  Bates  by  FSR  Spectroscopy."  FRK-116,  Departaent  of 
Oieplstry,,  Boston  College,  Cnestnut  HUl,  Massachusetts 
(October 


Ultee,  Casper,  Private  CoMMuiii  teat  ions . 

Vincent i,  ft.  C*.  , and  Kruger,  C.  M.  Introduction  to 
Phys  icail  Ga_s  nynawic.s , Kfw  York,  ftl ley  (isft^) . 

Leonard- Jones , J.  F;.  '"Processes  of  Adsorption  and 
Diffusion  on  Solid  Surfaces."  Trans.  Faraday  Society, 
No.  21,  p.  554  (1952).  “ 

Veast,  B.  C.  Mandiiiook  of  thewlstiy  and  Physics  (S5rd 
Edition),  t 1 e veTanT,~CKeai  »’c  al  Rubce  irTubl  j shi  ng  Coo  • 
pany  (19^2). 


Bardeen , 
Metallic 
(October 


J.  "The  tnage  and  Van  der  ftaal  Forces  at  a 
Surface."  Physical  Beviea,  Vol  SB,  p.  727 
194C).  


Nargenau,  H,,  and  Pollard,  ft,,  C.  "The  Forces  Between 
Keutral  Molecule*  and  Metallic  Surfaces."  Phys ica 1 
Ben  ew , Vol  oO,  p.  i2t  (July  1941). 

Prosen,  F,  >.  B.,  and  Sachs.  R.  C.  "The  Intrrnction 
Between  a 'loletulr  and  a Velal  Surface."  Physical 
Review , Vol  61,  p.  6S  (.lamia ry  1942). 

Mavroyann i 5 , C,  "Ihe  Intcr.oc?ion  of  Neutral  Molecules 
with  l»i»lrctric  Surfaces."  Molecular  Physics.  6,  p.  S9.t 
(196.5).  


Logan,  R.  M,,  ami  Fe^l  , . r.  "Classiral  Therrv  for  the 

Interaction  of  Cur-  with  Solid  Sur fines,’'  Journal 

of  Oie»^icaJ  Phj!?-,i^.«  Vol  4o,  *,o.  7 ji.  BoO  (.Uil>  . 

Pagni  , r,  I.,  and  Kivk,  J.  f.  "Diffusion  Throrv  for 
Adsorption  and  Desorption  of  f.as  .ttows  at  Surfaces. 

The  .lour  Mai  of  rherital  Plivsics,  Vol  50,  No.  5,  p.  1162 
TTebrmVfv  Tt!^) ‘ 

Hodak,  V.  T.,  and  Pagni,  P.  J.  "Ato«  Trapping  on  Sur- 
Fates."  The  Journ.il  of  rinrie.iS  Physics,  Vol  59,  No.  4, 
p.  20 n r^iguTt  i'5''.5r;  ' 


79 


Bibliography 

Chri nann  , , Schobtr,  0,  ^ F.rt  I , n . . and  N'cuttann » M. 

"Adsorption  of  llyilroprn  on  NicKcl  Single  Crystal  Suf' 
faces."  The  Journal  of  Choiaital  J’hysics,  VuJ  bO,  No.  11 
«.  4S2I  nunriTyrr, 

Gladstone,  S.,  I.aidler^  K.  J.,  and  tyvinn,  H Theory 
of  Rate  Process P'  . St*t  York,  McGraw-Hill  (l**^lin 

Van  dor  Ziel-  A.  Solid  State  Physical  Fleet  ronio . 
Eftflewood  Cl  iffs  , ?vev  Jersey,  T^rcnTTc^e'-TluTil  tT9TT) . 

SMith,  W,  V,  "The  Surface  Rcconbinat ion  of  H -Vtoas  and 
OH  Radicals."  Journal  of  Cheaical  Physics,  Vol  11, 
p.  110  (March  1 9 1 iT.  .v  . - 

Kursitjs,  S.  C.  "la*p  Reaction  Models  for  Analv^is  of 
Cheaical  f.ascrs."  TR-RK  "*R  5 ?1  i w i. t*s« 

iiirectoratp , ’J.s:.  Aray  Missile  Research,  I'ovelopaent 
and  rngincering  I ahoratory,  Redstone  ir^cnal,  Alabaaa 
(June  I'i-S}. 

Mi rsthff Idcr , J,  0,,  Curtis?i,  C.  r , and  Bird,  R.  B. 

Molrcul-r  llicoiv  of  Car.cs  and  iiuuids.  Nea  York, 

— 

Svehla,  R.  "I*t»natcd  Viscoritirs  and  Thcm.-il  Con- 
ductivities of  G.i'.rs  at  Mil'll  1 rwpe  I at  uro' . " NASA  TR- 
R-152,  Ni’tjonal  Ynonautii'-  and  ’ip.ict  Adwtn i s 1 1 at  ion 
(196?). 

Micklev,  M.  S.  , Ho*  • , R.  C.,  Sqmvrr'- , A.  1.  , and 
Str’.irt,  W.  r,  "Meat  Ma*'-  .ir.d  Mo-ii'ist u"  Tr.in-for  for 
Flnv  evn  a 11. it  Plate  with  Plowijj)'  or  t i <»n  . " 

NAC  Y • 1 N ■ Yr<’ .S  , Natien.dl  Av’vi'.ory  { Oi' -’i  1 1 1 ee  for  Aero* 

nai  t i V ( 1 F'  liA  ) . 

Pays,  K.  '1.  Cctw' v t ’ ve  Me.it  .ind  "a*.-  Tran‘fer.  .San 

Franc  I '■V  o , 'irCi  .o.'- IHIT 

'lyrrsoti,  A.  I . "Mcchat\  1 m-!’*  of  *^urf;(Ce  Keconh in.it  ion 
Itop  Step  r\;nitivn  Mow*,  o'l  Atovpji.  {‘'nv^en  (Aver  Noble 
Metal*-."  .'•ooin.il  of  ChcrpivHl  Plo^ics,  \ol  4?,  No.  9, 
p.  AJrii  (Hay-IWYT.'  ■■  


Appendix  A 


Sol  lit  ion  of  f ions  of  ^iotion 
for  the  To)  I is  ion  Model 


The  Kas'vall  collision  system,  as  shown  in  1 ijnire  A.l, 
is  analysed  as  .ollows. 


Mpurc  v.l.  r-as-h/tij  (olli'ioh  ''vsrer 


At  thr  iristiirl  the  collision  hCf'.jn*-,  the  ivn','-  I'f  action 

ttreoTM' 


n j’  ♦ « 1;^)  : • « n (A.l) 

»i  V k : * I V * 0 (A.  2) 

K r R 
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The  solution  of  Iquations  (A.l)  «nd  (A.!)  proceeds  as 
follows:  The  $olv.iion  for  Y and  2 Is  assumed  to  have  the 

fom 


t(t) 

• A Sin{«t) 

(A. 5) 

2(t) 

• b Stn(bt) 

(A.«) 

In  the  usual  manner,  the  assumed  solutions  are  differentia 
ated  for  the  first  ard  second  derivatives  and  the  results 
substituted  into  Equations  (A.l)  and  (A. 2).  The  result  can 
be  expressed  in  matrix  form  h)> 


2 

(■m  w*  ♦ k > • k 

• • 

A 

1 » « 

■» 

1 ♦ <k^  * k^)) 

A 

0 


(A.S) 


Thus  for  t non  trivial  solution  of  A and  A the  matrix  on 
the  left  must  he  sinpulrir.  Settinf:  the  drtcrninant  to  rero 


♦ k k "0 
R s 


(A.h) 


Solving  for  the  roofs  of  liquation  (A.(») 
f 


'1.? 


(m  k *0  k k 1 * ».  *«,  k 4m  r k k 


1/2 


(A.  7 ; 


12 


r 


i 

i 

1 


wherr  , rrfcrs  to  the  solutio;i  of  v.  with  thi'  ind  u, 

1 A I 

refers  to  the  ■iolvition  cf  vith  th*'  The  araj*  1 j t uJe  I 

ratio  e.^.’3tton  is  tS'en  by 


(K  - k )A  ■ It  II  - 0 
k k k 


which  gives  the  followins;  anpliiuJr  rat»"s 


(A.«) 


(A. 9) 


Thus  thr  g^r^^rll^  s<,jutipn  of  rqo.ltions  J\.l)  tiru!  (V,Z)  is 
g i ve,n  V v 


Mtl  ' •.j'^»n<-jt  * • A,s»nf,,t  ♦ ,„i  (A.lhi 

7 ( t ) • h , 0 f n : , T • . , ^ • H , s i n 1 , , I * I ^ ^ ( A . 3 n 

I I I A 4^  a. 

Iqiiiitirji'-  I ■'i . ! •/ ' .«:u'  .(’.11'  li.tvo  4 ii  MiU  j't’Ju^cn  t ursl  ru'-vti 


whit'h  .Tif  t’'».  init  i>:!  i ( iiJ  , t i hh  s 


Y((>)  - .•((■<)  rSui' 


(A . 3 :) 


vhcrc  C I*'  th'.'  i (li  t ! .1 1 ohasA  •nu'le  .iih?  i is  thr  I i t uile  of 


surf.n  o .ilor  use  I 11,1 1 son. 


t(0)  - 

COMi 

(A, 

t(0)  • 

' »' 

CA. 

wtiera  c ia  th«  vnerfy 

1 / 

/rar jnrpT  \ 2 

^ * r[h 

•/v' ) 

(A. 

•s  indltatriii  In  rquatlon  (1 

of  the  ■ain  text. 

let  the  anplitudir  ratios  frt>m  Tq^uation  (>\.9)  be 

denT?r*d  as  follows 

*1 

(A. 

^2 

„ ^2 
IT, 

(A. 

fhcn  the  r'snlt  of  vising  the  initial  ? 'Ttndi  t ion*,  yields  the 
solut ions 


Y(t)  - (A  jro<V  I ) i n«*'j  t * ( A ^Sini  ^ )cost,j  t 

♦ (AjCosi  j)Sinwt2l  ♦ (A^RIrti  (A. 


Z(t)  • AjCOs^j ) SiBtrf  1 1 ♦ (A^Sinf  |)co5tu^  ;i 
1 ' 

♦ ^ |(Aj,co»^j)Slnv2*  * (A^Sinl  2)cwistoj  :| 
2 

?(t»  - j )<'o^«''|t  - y,  CAjSini  j >Sinwjt 


whcrr 


AjSlnt^j  - “(0) 


A^Sini:^  * rtf') 


- (-1^ 


A,co«.^  • (:~) 


* 1 
* 1 

^■2 

T""'"  '•  T 
r ■ r 

' I 


1 

T 1 

' 2 ' 1 

^*Y(o)  - :co) 

■ T 

- r: 


2(0)  ‘ j-M'(0 

r ■ r 
'2  ^ 1 


The  results  of  coBj»*itat  ions  wruie  fron  i qiintions 
(A.  19),  and  (A.  70)  for  a t e"'per  a t ui  e of  ^00  H and  an 
non-dimensional  velocity 


(A.  19) 


(A. 20) 


(A.  21) 

(A. 22) 

(A. 23) 

(A. 24) 

(A.  1«)  , 
inconi  nj* 


' 0.5  (A.25J 

for  fluortnv  on  copper  Is  (iven  in  Figure  A.  2.  Exeiiinatlon 
of  Figure  A. 2 shows  that  for  a nondiocnsional  initial  v«- 
ocitjr  of  0.5  the  probability  of  being  trapped  (represented 
by  *'i^**2  Figure  A.  2)  Is  50  percent. 

CoMputations  aiade  for  a range  of  initial  velocities 
yields  Figure  A. 5 (note  in  Figure  A. 5 that  the  protibility 
corre'iponding  tn  an  initial  non  dineesjonal  velocity  of  0.5 
is  O.S). 

A faoily  of  curves  si«llai  to  the  curve  In  Figure  A. .5, 
for  a range  of  teaprratures  fro*  0 K to  500  t anJ  Figure  4 
froo  the  oain  tr*t,  coobine  to  yield  the  curve  for  the 
initial  st^cVing  coefficient  such  a the  broken  line  in 
Figure  9 of  the  naln  ,e»t  for  the  stteVing  coefficient  of 
fluorine  on  copper. 


Hu 


0,0  0.:  o.f  o.«  1.0  1.4  \.n  j.f, 


IHITI/k.  MilZ  AT  or  CTALlSlOri  f (R;\PUiS  X |t) 

riftuTi  A,'.  tnltfal  AmcJo  at  On>>rl  of  COili^ton  v.<, 
Afpn-iiut  IfiirRS' 


Mp.urt'  A.V  I’rohalu  1 1 1 V of  Atop  He  ini;  Tr.'U’pcJ 


Appendix 


Pt  f r»nce  > qu*t  toiffs 

All  functions  o<  x,  and  f y»trt  treated  as  the  follovinR 
fictitious  function  f(v,t)  usinic  the  cell  nooenclature  shown 
in  Figure  H . 1 . 


♦ ’*k' 

k’  k*’'  k 

< *'  < '-k‘  o V'  k 

k k* 

i « ’ ■ k c'  (I  S’*  (k< 

HfM  I 

I i purr  H,!.  Tel  N’onrn*.  Imt  urr 
lif'l  ami  «ccond  Jrriv,'itvrs  aj'pt  o*  rd  h> 


a f ( X , i , M ( s , , ) 
>\ ■ 


■ ‘-T  * ■' 

tiheif  tfie  I*  .iM  defmrd  in  th«’  fcik'ittini 


FIRST  DERIVATIVE  A OPERATOR 


Sl|i 


Af(x,4»^)  - Bf(x,i^g)  - Cf(x,i{>^) 
D 


(B.3) 


where : 

for  all  cells  except  the  wall  cell  and  the  last  cell 
A = 

B = k’^ 


c « 

k^  - k’^ 

D = 

k^k'  + kk’^ 

' i|/'  + k' 

- k 

^c- 

£or 

the  wall  cell 

A = 

(k’  + k”)^ 

B = 

k’^ 

C = 

(k'  + k”)^  - k’^ 

D = 

k’^k"  + k'k"^ 

*A  = 

il;’  + k' 

!{,’  + (k'  + k”) 

for  the  last  cell 
A = k 

B = (k  + k’)^ 
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C « - (k  + k’)^ 

D « k^k’  + kk'^ 

- k 

ijjg  » ij;'  - (k  + k') 


SECOND  DERIVATIVE  A OPERATOR 


Alp  ^ 


^=ip' 


jAf(x,ipj^)  + + Cf(x,4)^) 

- ; B 


(B.4) 


where ; 

for  all  cells  except  the  v/all  cell  and  the  last  cell 
A = k’ 

B = k 

C = -(k  + k’) 

D = k'k^  + k'^k 
= 4)’  - k 
4)g  = 4^'  + k' 

= r 


for  the  wall  cell 
A = (k'  + k") 

B = -k' 

C = -k" 

D = k’^Ck’  + k”)  - k’ (k'  + k”)^ 
4^^  = 4;'  + k' 


90 


• t’  • (k*  ♦ k") 


lor  l«»t  cell 
A - (k  • k’) 
li  - -k 
C - k’ 

n • (k  M k’)i^  - (k  ♦ k')*k 

'♦a 

- (k  ♦ kn 

♦c  • ♦' 

DIFFrRIVTr  EQUATIONS 

ll»lnE  the  A opcrutors  defined  above.  Fquatlonft  (AS), 
(46),  and  (4?)  vere  approx iiaated  for  the  computational 
scheme  as  follows. 


ui«*  X)  - u{.)  . X,  j.  -i  g 

• 4 ('■‘"'I'v  ■ * ™",7  ! 


r.(x.x.i  . .j(«i  . ;.x  U 


A * ^1  I 

u f'i'")  fr  ■ ’j'”’ 
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T(x»A*) 


fcpu 


II 


Hi^  pressurr  wns  fitted  «'ith  an  analytic  curve,  thus 
was  taken  frosi  the  derivative  Of  the  analytic  curve. 


(i.7) 


dP 

r* 
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Appendix  C 


Drivat  i'on  of  Aomtdary  Condi  lient  for 
£ Concont ration* 

TIkoro  are  t*fo  way*  of  approacMng  the  derivation  of  the 
concentration  boundary  condition*.  One  way  is  purely 
analytic  and  will  he  diseuxsed  first;  another  way  is  to 
exaaiine  the  analytic  derivation  is  retroispect,  thl*  will  he 
done  second. 

ANAtYTIC 

Since  the  no  slip  boundary  condition  for  velocity  is 
used,  it  would  not  seen  a had  assumption  to  consider  the 
region  from  the  center  of  CELL  2 to  the  wall  as  hein^ 


Figure  C.l.  Conservation  Fquations  for  Stagnant  Cell 
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essentiallx  statnant.  Further  the  horiiontal  concentra< 
tlofi  gradients  arc  as^uated  to  be  TuegHgJble  toaipared  to  the 
vertical,  these  atsuaiptiont  lead  to  an  analysis  as  shown  in 
In  Figure  C, 1 . 

The  conservation  of  species  equation  in  steady  state 
is  then  written 


37  ® 


iv  Fick's  law 


0J.V, 


dOr 


Thus  the  governing  equation  Is 


j U 

Ty^fTf} 


Intcgratlun  yields 


dop 

Pp  • A 


(C.l) 


(C.2) 


CC.^) 


(C.4) 


where  ^ is  a constant.  In  partirulai,  at  the  catalytic  wall, 
y • I.,  the  steady  stale  flux  oust  be  equal  to  the  rale  of 
recombination  and  attachaient: 
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DiffuKtv*  flux  ^ 2tN 


(C.5) 


wlitre  the  2 accounts  for  thr  replacoairnt  of  the  uall  aio* 
4«tacht>d  by  tha  recooblnat ion.  Tbis  it  nacaisary,  to 
Maintain  the  o^ui  1 ib*^iun  wall  concentration  of  F atoaii. 
Thu* 


2y« 


1/2  M dnp(b) 
* * — 


(C.b) 


And  thus  A nuvt  be 


A • - Zyn, 


(C.7) 


We  n^'w  have  the  differential  equation 


dn. 


F “Uy” 


(C.I) 


Assunii'iK  P^,  is  corstant  and  inCeigt  at  ini;  yielvis 


n(  •) 


tnj.(i.)  / \i;: 

‘-TV“  \Tr^}  ^ • > 


(C.9) 


Since  at  y • 0 the  concentration  of  F mist  be  (n,  !» 

r CLLlrf  4. 

mist  equal  j-  Thus  we  hnve 


1 


r ♦ 1 (c.io) 


n|.-(y) 

W^CELL  2 


^ np(L) 

^PcEtt  2 


Finalljr  sine* 


»p(y) 

^F^CELL  2 


f I") 

^"f^CELL  2 


•nd  'ealUlnt  that 


^"r^lL  2 


1 ♦ 


n^*F 


lY, 


dn_  (I.)  V I / ? 


(C.Il) 


(C.12) 


(C.IS) 


sSnet'  fach  rcconUin^^t  Ion  ro^ults  in  the  cmtlon  of  »u  F^, 
»ol«co)r.  Then  in  a vinnnr  Nnnnrr  to  that  used  in  gottinf; 
Ffjuat  5 or  (C.  i) , 


dn 


f . 


’.■j  -sr  ■ 


• "’(r'irr* 


((.141 


• nd  intfR'ir.it  ion  yiolds 


"i:  (y)  “ 
2 


If".  <”(3%) 


1/2 


)•  ♦ B 


(C.IS) 


9fr 


Since,  «|«ln, 

^"F|*criL  r 


np  (y*0)  I#  (Hj. 
Thus  knowin;;  <n 


>CELL  2*  •'*** 

z 

F^WAU  Equation  (C.i2^ 


mall 


^*pVall 


^"F^^CBLL  2 


(C.16) 


WTMOSrECT 

Sti'/riint  fioM  Equation  (C.12) 


^"f^csli  2 


and  roarianginR  it  slir'.tly  yield;; 


(C.i7) 


[^"r^Al.l 


(n») 


F'etU.  21 


(C.II) 


Equation  (C.H)  the  identical  statrnrnt  to  as»  <Ming  that 
the  slope  of  the  concentration  at  the  wall  is  li.iear  and 
k'Pf rox inat  ing  the  derivative  by 


doj  ("r^KAi.L  ■ 2 

^ ^ 


(C.l?)) 


Thus  in  retrosfei't  one  could  have  derived  the*  boundary 
conditions  for  r and  by  ncrely  ns  uming  a linear  profile 

4 

at  the  waM  . 
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Appendix  D 


Teaperetuye  egd  Pressure 
Dependence  of  the  Fluorine 
Veil  Reconbinet ion  Coefficient 
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Hi  !^!io<c  i at  Ion 
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